arXiv:hep-ph/0211068vl 5 Nov 2002 


OSU-HEP-02-14 

UMD-PP-03-026 


Lepton Flavor Violation and the Origin of 
the Seesaw Mechanism 


K.S. Babu^, B. Dutta^ and R.N. Mohapatra^ 

^Department of Physics, Oklahoma State University, Stillwater, OK 74078, USA 
‘^Department of Physics, University of Regina, SK SfS 0A2, Canada 
^Department of Physics, University of Maryland, College Park, MD 20742, USA 


Abstract 


The right-handed neutrino mass matrix that is central to the under¬ 
standing of small neutrino masses via the seesaw mechanism can arise either 
(i) from renormalizable operators or (ii) from nonrenormalizable or super- 
renormalizable operators, depending on the symmetries and the Higgs con¬ 
tent of the theory beyond the Standard Model. In this paper, we study lepton 
flavor violating (LFV) effects in the first class of seesaw models wherein the 
vr Majorana masses arise from renormalizable Yukawa couplings involving a 
B — L = 2 Higgs field. We present detailed predictions for r —> /r -|- 7 and 
p, ^ e + j branching ratios in these models taking the current neutrino oscil¬ 
lation data into account. Focusing on minimal supergravity models, we find 
that for a large range of MSSM parameters suggested by the relic abundance 
of neutralino dark matter and that is consistent with Higgs boson mass and 
other constraints, these radiative decays are in the range accessible to planned 
experiments. We compare these predictions with lepton flavor violation in the 
second class of models arising entirely from the Dirac Yukawa couplings. We 
study the dependence of the ratio r = B{p e + 'y)jB{T —> p + 'y) on 
the MSSM parameters and show that measurement of r can provide crucial 
insight into the origin of the seesaw mechanism. 
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I. INTRODUCTION 


With the evidence for neutrino masses and mixings from solar and atmospheric neu¬ 
trino data getting more and more firm, the nature of new physics that could explain the 
observations is under a great deal of scrutiny |]l[. There are several issues that need to be 
understood, notably (a) the smallness of neutrino masses, and (b) the origin of the large 
atmospheric neutrino as well as the solar neutrino mixing angles, the latter being favored 
0 by the combined solar neutrino results 0 including the recent SNO neutral current data 

Seesaw mechanism [Q provides one of the simplest ways to understand the small neutrino 
masses. It assumes the existence of a set of three right-handed neutrinos {vr) that have 
masses at the scale vb-l corresponding to a new (local) B—L symmetry of weak interactions. 
Atmospheric neutrino observations suggest that the scale vb-l is much lower than the Planck 
scale, leading to a new threshold inferred solely from experiments. At low energies, these 
heavy right-handed neutrinos induce operators of the form /vb-l-, where W is the 

Dirac Yukawa coupling matrix connecting the pb with the left-handed neutrinos and 
the Standard Model (SM) Higgs doublet Hu- L denotes the lepton doublet which contains 
pl- After electroweak symmetry breaking, these operators induce masses for the light pl 
of order Y^^v'^^/vb-l- Since vb-l is much larger than the weak scale v^k, the induced pl 
masses will be extremely small compared to all other masses of the SM fermions. Present 
neutrino oscillation data suggests a scale vb-l 1012 _ iQis GeV. 

The precise seesaw formula for the light neutrino mass matrix with three generations is 
given by 


M, = -MIM],^Md , ( 1 ) 

where Mq is the Dirac neutrino mass matrix and Mr is the pr Majorana matrix. The 
neutrino mixing angles in such schemes would arise as a joint effect from two sources: (i) 
mixings among the right-handed neutrinos present in Mr and (ii) mixings among different 
generations present in the Dirac mass matrix Mr. The (physical) neutrino oscillation angles 
will also receive contributions from mixings among the charged leptons. 


A. Neutrino mixings and lepton flavor violation 

The neutrino flavor mixings induced by the seesaw mechanism (as needed by oscillation 
data) can lead to lepton flavor violating (LFV) effects. However, within the SM extended 
minimally to accommodate the seesaw mechanism, such effects are extremely small in any 
process other than neutrino oscillation itself, owing to a power suppression factor {I/vr-lY, 
as required by the decoupling theorem. This situation is drastically different if there is low 
energy supersymmetry. The power suppression then becomes {I/Msusy)"^ which is very 
much weaker and can lead to observable LFV effects at low energies, as noted in a number 
of papers The main difference is that with low energy SUSY, lepton flavor violation 

can be induced in the slepton sector, which can then be transferred to the leptons suppressed 
only by a factor (I/Msusy)"^- The way this comes about is as follows. Consider the minimal 
N = 1 supergravity (mSUGRA) models 0. At the fundamental scale Mpi where SUSY 


2 


breaking is communicated to the SM sector, all the soft SUSY breaking scalar masses are 
taken to be universal, and the trilinear A-matrices are taken to be proportional to the 
respective Yukawa matrices. Thus, at the Planck scale, there is no flavor violation anywhere 
except in the Yukawa couplings and SUSY-preserving mass terms in the superpotential. 
Now, the right-handed neutrinos have masses of order vb-l which is much lower than Mpi, 
as inferred from atmospheric neutrino data. In the momentum regime vb-l < At < Mpi 
where the up helds are active, the soft masses of the sleptons will feel the effects of LFV 
in the neutrino Yukawa sector through the renormalization group evolution. At the scale 
vb-Li the slepton mass matrix is no longer universal, and this non-universality will remain 
down to the weak scale. This LFV in the slepton sector is subsequently transferred to the 
leptons through one-loop diagrams involving the exchange of gauginos. 

It should be noted that in the absence of neutrino masses, there is only one leptonic 
Yukawa matrix, for the charged leptons, which can be diagonalized at Mpi. will 
remain diagonal to the weak scale, and would not induce any flavor violation in the slepton 
sector in mSUGRA models. The experimental evidence for neutrino oscillation is thus a 
strong indicator that there might very well be lepton flavor violation, assuming the validity 
of low energy SUSY. Searches for LFV processes such as r —/r + 7 and/or n ^ e + 'j can 
therefore be an important source of information on the up mixings in Mp and/or family 
mixings in Mp. (There are some exceptional circumstances, which we point out later in this 
section.) 


B. Majorana LFV versus Dirac LFV 

It is possible to classify the seesaw models into two classes depending on the way the vp 
gets its Majorana mass Mp. In one case, Mp arises from renormalizable Yukawa couplings 
involving a B — L = 2 Higgs field A through the superpotential terms 

W D /z/V"A + Y^vMHu . (2) 

Here / is the Majorana Yukawa coupling matrix, Y^ is the Dirac Yukawa matrix and we use 
the standard notation appropriate for supersymmetry with being the conjugate of up. In 
the second class of models, the matrix Mp is either put in “by hand” as a bare mass term 
in the Lagrangian (such terms are super-renormalizable), or arises from non-renormalizable 
operators involving a B — L = 1 Higgs boson via the couplings /^■ models of 

the second class, the flavor violating slepton mixings can arise (assuming mSUGRA) only 
from the Dirac Yukawa couplings of the neutrinos, Y^. On the other hand, in the hrst class, 
which is the main focus of this paper, LFV can arise from both the Dirac Yukawa coupling 
Yiy as well as the Majorana Yukawa coupling /. 

For the purpose of the present paper, we will call the two lepton flavor violation alter¬ 
natives as LFV in the Majorana case and LFV in the Dirac case. 

The simplest example of models with renormalizable Majorana Yukawa couplings is 
provided by the minimal (SUSY) extension of the SM with a gauged B — L symmetry. 
Bare masses for the vp are then forbidden hj B — L gauge invariance. If the Higgs sector 
contains a B — L = 2 field A, the Yukawa coupling /z/^z/'^A as in Eq. (2) will be allowed 
which will induce Mp of order vb-l = (A). (If instead of A, a R — L = 1 Higgs held 
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is used, Mr can arise from a non-renormalizable coupling /^■ This belongs to the 

second class.) There are interesting extensions of this class of models with renormalizable 
Majorana couplings, such as SO(IO) with a 126 Higgs held or left-right symmetric 
model with a triplet Higgs held. In these models the Majorana Yukawa couplings may play 
a dominant role in generating the required neutrino mixings. The role of the Dirac Yukawa 
coupling Yy for LFV may be subleading and may even be negligible. To be concrete, we 
shall assume in our analysis that the matrix is diagonal in a basis where the charged 
lepton Yukawa matrix Yi is also diagonal. For numerical purposes we shall further assume 
that Yjy and are proportional, a situation realized naturally if the gauge symmetry is left- 
right symmetric. The neutrino mixings will all arise from the Majorana coupling matrix /, 
which will be determined (upto an overall scale factor) from neutrino oscillation data. The 
proportionality assumption Yy (x. Yi is not crucial to our main conclusions on LFV ehects 
related to Mr, but it greatly simplihes the presentation of our results. This assumption may 
be viewed as analogous to the assumption often made in the case of Dirac LFV models that 
the super-renormalizable Majorana mass matrix Mr is proportional to an identity matrix. 


C. Summary of results 


In an earlier paper we have presented results of our preliminary investigations on 
LFV involving the Majorana Yukawa matrix /, assuming proportionality of Yy and 1^. (We 
called this relation up-down unihcation.) There we studied a general class of SUSY breaking 
models and did not adhere to mSUGRA. The analysis of Ref. ii was carried out in the 
context of left-right symmetric gauge theories with an enhanced gauge structure. In contrast 
to Ref. |T^, in this paper we analyze LFV effects induced by the / matrix within mSUGRA 


models. This means that the only source of lepton flavor violation is in the neutrino Yukawa 
sector. Signihcant LFV effects are found even in this minimal scenario. As explained in 
Sec. H.G, a novel way of htting the low energy neutrino oscillation data has the effect 
of enhancing the LFV in this minimal scheme. We also use a simplihed gauge structure 
involving only an extra B — L symmetry. Even this symmetry is not essential, the crucial 
ingredient is the Yukawa coupling matrix / of Eq. (2). Thus the models studied here are 
more minimal compared to the left-right models and are likely to be the low energy limits 
of a wider class of unihed models. 

Our main result is that LFV effects associated with the Majorana Yukawa coupling / are 
in the range accessible to forthcoming experiments for a large range of MSSM parameters 
consistent with constraints from dark matter, the Higgs boson mass and direct experimental 
search for supersymmetry. We also point out that the detailed dependence of the branching 
ratios in processes such as r —p -|- 7 and /i —> e -|- 7 on the MSSM parameters is different 
in this class of models compared to the case of flavor violation associated with the Dirac 
Yukawa coupling Yy. With some information of the MSSM parameters, measurements of 
these radiative decays will enable one to pin down the nature of interactions that give 
rise to the right-handed neutrino masses. Obviously, this will have a profound impact on 
our understanding of the neutrino sector and can shed light on the origin of the seesaw 
mechanism. 

It may be worth pointing out that there are other experimental predictions that may help 
to distinguish between these two possibilities. For instance, two of us have recently pointed 


4 
















out that the Dirac LFV models can lead to an observable signal in the baryon number non¬ 
conserving process neutron-anti-neutron oscillation |]^] whereas in the Majorana LFV case, 
the N-N oscillation is unobservable. 


D. One caveat 

There is one caveat in the foregoing discussions that should be mentioned. It applies 
equally well to the Dirac and Majorana LFV alternatives. In a special class of seesaw schemes 
it is possible that no lepton flavor violation shows up in the next round of experiments 
without contradicting the current neutrino oscillation data even in the presence of low energy 
SUSY. While we feel that such situations would be non-generic, it should be born in mind, 
nevertheless. In the Dirac LFV alternative it occurs as follows. It is conceivable that all the 
flavor mixings arise from the mass matrix Mji and that Y^, is diagonal. Flavor violation in 
the bare mass parameter (or non-renormalizable term) Mji is not transmitted to the slepton 
sector via the RGE, implying highly suppressed LFV effects. One can easily £t all the 
neutrino oscillation data via the mixings present in It is also possible that the scale 
Mu is much below 10^^ GeV, in which case the elements of will be much smaller than unity 
so as to fit the neutrino data, again suppressing LFV effects. Similarly, in the Majorana 
LFV alternative, it might so happen that the elements of / are very small, compensated by 
a higher value of vb-l so that Mb is unchanged, or that the elements of V, are small. Such a 
scenario will suppress lepton flavor violation. A null result in the rare lepton decay searches 
would point to a situation in which either of these two situations could be operative and we 
will have no way to tell which operators are inducing the vb masses. On the other hand, 
if there is a positive signal in the rare decays, our results will help in probing the origin of 
Mb- It appears to us natural to have at least the third family Yukawa coupling in Yy of 
order one (analogous to the top quark Yukawa coupling being of order one) and similarly at 
least one element of / (/ss) to be of order one. Observable LFV effects follow in this case, 
which is what we analyze. 


E. Outline of the paper 

We have organized this paper as follows. In Section II, we present the basic model and 
the relevant renormalization group equations (RGE) that will be used in getting the slepton 
mixings and thus the lepton flavor violating rates for the case of Majorana LFV. Here we 
also present analytic approximations to the LFV effects by integrating the RGE. These 
expressions, while only approximate, are useful in gauging the signihcance of LFV in the 
Majorana case as well as in comparing it with the case of Dirac LFV. In Section II.A we 
present the relevant expressions for obtaining the / matrix by inverting the seesaw formula. 
In II.B this information is used to compare analytically the predictions for r —> /r -|- 7 and 
p —> e -|- 7 branching ratios in the two cases. In Section II.G we present a convenient way 
to parametrize the Majorana Yukawa matrix / that can be compared directly against low 
energy neutrino oscillation data on the one hand and LFV effects on the other. Section III 
is devoted to our exact numerical results for LFV effects in the case of Majorana as well as 
Dirac Yukawa couplings. In UFA we present our numerical hts to the neutrino parameters. 
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III.B has our results for the branching ratios r —/i + 7 and /r —e + 7 as well as the ratio 
of the two branching ratios. Section IV has our conclusions. 

II. SLEPTON MIXING AND LEPTON FLAVOR VIOLATION FROM 
MAJORANA YUKAWA COUPLINGS 

In order to illustrate in detail the source of LFV for the case of Majo- 
rana Yukawa couplings generating the vr mass matrix (the Majorana LFV alter¬ 
native), we work with a simple extension of the minimal supersymmetric stan¬ 
dard model that accommodates the seesaw mechanism. It is based on the 
weak gauge group SU{2)l x U{l)i^^ x U{1)b-l with fermions assigned as fol¬ 
lows: Q( 2 , 0 ,+|); L( 2 , 0 ,-l); «'=(!,-i,-i); +i,-|); e'=(l,- 1 -|,+ 1 ); -M). 

This is the simplest extension of the weak interaction sector that predicts the 
existence of the vr from gauge anomaly cancellations conditions and thus to 
small neutrino masses via the seesaw mechanism. The Higgs superfields are 
iL„(2,-|-i, 0); iLrf(2, —1,0); A(l,-|- 1 , — 2 ), A(l, — 1 ,-|- 2 ). The superpotential for this the¬ 
ory is: 

W = YuU^QHu + + Y.e^LHd + (3) 

+ /z/V'^A + + 5(AA - M^) . 

Here all the Yukawa couplings are 3x3 matrices; S' is a singlet held introduced to help 
facilitate the breaking of U{1)b-l x U ( 1 )/ 3 jj down to U (l)y in the SUSY limit. One identihes 
the SM hypercharge as Y/2 = + {B — L)/2. This superpotential gives equal vacuum 

expectation values (VEVs) to (A) = (^A^ = vb-l- The Yukawa coupling matrix / will 
induce large mass for the held, leading to the seesaw mechanism. This model should be 
contrasted with the commonly discussed case in the literature where the heavy vr masses 
are either put in by hand or arise from a non-renormalizable term in the superpotential. 

In order to study havor violation, we will make the simplifying assumption that both 
Yu and are diagonal at the Planck scale; therefore all neutrino mixings arise from the 
general havor structure of the Majorana Yukawa couplings fij. Then from the approximate 
bi-maximal pattern for neutrino mixings deduced from oscillation data, we can obtain 
for diherent patterns of neutrino masses. Since the model is supersymmetric, the havor 
violation present in fij will induce processes such as r —/i -|- 7 and /i —> e - 1 - 7 decays. We 
investigate these predictions using simple analytic approximations in this section. 

We need to write down the RGE for the Yukawa couplings and slepton masses for the 
class of models under consideration. Denote the left-handed slepton mass squared matrix as 
and the trilinear A terms in the slepton sector as A^, Af (in an obvious notation). 
The relevant RGE for Y^, and mf are: 

ctY Y 3 

'W “ i^l^''(3U7 + u7) + 3yX + y;y, + 4/'/-3!,= -!,J--9yj (4) 

W iYy~i " T 

" ley 2’"b)y/y< + ("»i + 2™k)yju + 2Y!mlY, + YiY,ml (5) 

+2y>;,y„ + yJUmi + 2431 , + 2AIA, - 6j|M| - 39|_iAf|_J 
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The RGE for the relevant A parameters are: 


^ = -^yMTriWlY, + Y/n) + hY}Yi + YJY, ( 6 ) 

- -si- IsLlI + Y,lTr(6A,Y^ + 2A,Y}) 

+ 4y/ + 2Y^A^ + Qg^M2 + 2g\MR + 35'g_^Ms-L]] 

^ + yjy;) + 5 yjK„ + y/y, (7) 

+ 4 /V - 39I - sj - ?j|_ J + K|rr( 6 yl„yJ + 27 l„yJ) 

+ 41 ^^ 24 ,^ + Sf^Af + 2Y^ Ai + Qg^M2 + 2g\MR + 

We are interested in LFV arising primarily from the matrix /. Adopting mSUGRA, we have 
at the Planck scale, = A^Yy, Af = Aof. Furthermore, at Mpi, we have 

universality of all scalar masses (denoted by mo) and a common gaugino mass (denoted as 
^^ 1 / 2 )- 

Under the set of assumption specihed above, it is clear from Fqs. (4)-(7) that the 
off-diagonal elements of / will not induce slepton flavor violation to the lowest order in a 
perturbative expansion, since / does not enter directly in the RGF for m| or However, 
the one-loop solution for m| involves YjYy, which receives a contribution from /!/ (see Fq. 
(4)). Therefore in second order in the parameter we do have flavor violating 

effects in m| and A^. It is important to note that this is not a two-loop RGF effect, rather 
that it is a one-loop RGF improved effect, as indicated by the presence of the [in (]\^^)]^ 
term. 

We can estimate the strength of the slepton flavor violation by integrating the RGE 
analytically. We hnd from Eq. (4)-(7) that 

# j) ^ zYA^lYXf'f + ffYXX («" ( 8 ) 


itj) ^ ^{AXXXf + /'/yjyji.i (<"-^) 


The branching ratio for ij —> + 7 can be written approximately as 

B{ij ^ ij + 'y) ^ ^aemjcigl + gjf / Am| y 2 

B{ij ^ ii-\- Uj Pi) 327rm‘^iGp \ mh ) 


( 10 ) 


Here mh is the average slepton mass at the weak scale, which may be quite different from 
mp. Eq. (10) is obtained as an approximation to the chargio and neutralino exchange 
diagram proportional to LFV in Eq. ( 8 ) with the enhancement factor of tan (3 as indicated. 
Cl in Eq. (10) is an 0{1) coefficient. LFV arising from Eq. (9) can be comparable, but 
it will not change our approximate results by much. We have found that Eq. (10) is a 
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reasonable approximation to the exact numerical results to within a factor of few, provided 
that tan/3 is not too large and that there is no large hierarchy between the SUSY breaking 
mass parameters (p, ttIsi, 'rnii 2 )- 

Within mSUGRA, mi /2 niust be larger than about 250 GeV, to be consistent with the 
Higgs boson mass constraint (viz., mh > 114 GeV) and the b ^ s'j constraint. Since we 
also want to identify the lightest neutralino (LSP) as the dark matter in the universe, the 
relic abundance puts another constraint that the second lightest SUSY particle (NLSP), the 
right-handed stau in our case, should be nearly degenerate with the LSP to within about 
20-30 GeV. This condition can be satished if mi /2 — 4.4mo is obeyed, which we shall adopt. 
Take for example a relatively light SUSY spectrum with mi /2 = 250 GeV and tan/3 = 10. 
The charged Wino mass is then approximately 200 GeV, and the Bino mass is about 100 
GeV. Ghoosing mo = 55 GeV and Hq = 0 leads to a mass of about 190 GeV for the 
left-handed sleptons and a mass of about 130 GeV for the right-handed sleptons, with the 
right-handed stau being somewhat lighter around 115 GeV. Go-annihilation of LSP will be 
efficient with such a spectrum. 

Next we need to choose the neutrino parameters. With V, oc a good £t to the current 
neutrino oscillation data can be obtained with the choice = 0.66 and (/V )23 — 0.28 

(see detailed numerical fits later). Furthermore, in{Mpi/M b-l) — 9 is quite plausible. With 
this choice, we estimate B{t —/U -|- 7 ) ~ 2 x 10“® for tan/3 = 10 from Eq. (10). As Aq 
is increased from 0 to 300 GeV, this branching ratio increases to about 2 x 10“®. Thus 
we see that for interesting ranges of parameters, the process r ^ fi + 'j can be within 
experimental reach. More careful examination of the decay is thus warranted. Similarly, 
the decay /U — e -+- 7 can have a branching ratio of order 10 “^^ — 10 “® for reasonable ranges 
of model parameters. Some of this parameter space is already excluded by data, there is a 
good chance of being able to measure this decay in the planned experiments. 

It is instructive to compare the branching ratios obtained here in the case of Majorana 
LFV to that expected in the case of Dirac LFV. In the case of Dirac LFV, assuming that 
the Up Majorana mass matrix is proportional to the identity matrix, we can write down the 
analog of Eq. ( 8 ) as 

^ j) ~ -^(3mg + Al){Y^Y^)ij ■ (H) 

Gomparing Eq. (11) with Eq. ( 8 ), we find that numerically the two effects can be compa¬ 
rable. Which effect is more dominant will depend on the overall strengths of the coupling 
matrices in the two cases. This will be discussed further in Sec. III. 


A. Flavor structure of the Major ana-Yukawa couplings 

We can invert the seesaw formula (Eq. (1)) and obtain the matrix elements of the 
Majorana Yukawa coupling matrix /: 

/ = . ( 12 ) 

vb-l 

We also have the relation for the inverse of the light neutrino mass matrix in terms of 
the light neutrino mass eigenvalues (mi, m 2 , m 3 ) and the leptonic mixing matrix U as 



(13) 


M 


-1 

u 


( ^ 



0 

0 


0 

m 2 ^ 
0 


0 

0 

mg^ 


U^. 


We ignore CP violation and nse the near bi-maximal approximation for U given by 


U ~ 


^ c 

s+ce 

x /2 

s—ce 

\ V2 


s e ^ 

c—se 1 

V2 P 2 

—c—se 1 

p 2 V2 / 


(14) 


. s is the 

solar nentrino mixing angle for which the allowed range after the SNO data is [0 0.71 < 
sin^26*0 < 0.94 at three a conhdence level. Mnltiplying all these together, we have 


Here Ue 3 ~ e < 0.16 from the CHOOZ and Palo-Verde reactor experiments Id 




(16) 


where we have assnmed Md to be diagonal with its entries denoted by mD,i- The dehnition 
/i“d = has been nsed. We have then 


h-ii 

/^12 

/^13 

^^22 

1^23 

/^33 


£2 

- 1 - 1 - 

mi m 2 m 3 

c(s-f-ce) s(c — se) 

^ ^ ^ + 


'/2mi v^m2 V^Tn3 

c{s — ce) s(c - 1 - se) 

y/2mi \/2m2 V^r«3 

(s -1- ce)^ (c — se)^ 1 
+ ^ + 


+ 


2 mi 


2 m 2 


21713 


(^2 _ ^ 2 ^ 2 ) (^2 _ ^ 2 ^ 2 ) ^ 1 


2 mi 


27719 


2m3 


(s — ce)^ (c-hse)^ 1 
+ - - + 


2 mi 


27719 


27713 


(16) 


Using these expressions, for different scenarios of nentrino masses we can determine the 
valnes of fij which enter into the flavor changing slepton masses. We will focns primarily 
on the hierarchical nentrino mass schemes where mi -C m 2 -C m 3 . 


B. Comparison between Majorana LFV and Dirac LFV alternatives 

Let ns hrst discnss the case where the right-handed nentrino masses arise from a renor- 
malizable Majorana Ynkawa conpling /. As in the previons section, let ns assnme that the 
Dirac nentrino mass matrix Mu is diagonal in the same basis where Mi is diagonal. One 
interesting pattern is provided by the normal hierarchy mi -C m 2 -C m 3 . In this case, it 
is plansible, bnt not mandatory, that in Eq. (16), the terms with mi in the denominator 
dominate /i“^ with all other terms being negligible. In this approximation, we obtain for 
the process r —/i -|- 7 process an amplitnde proportional to 
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(17) 


(Am^)23 oc 


mD,2m%3s'^ 

Amlwlkvl_L 


Similarly, the /x —> e + 7 process has an amplitude given by 


(Am^)i 2 oc 


2S^c 


(18) 


We have set e = 0 for simplicity here. From Eq. (17)-(18) it might appear that the branching 
ratios can increase without limit as mi is decreased. However, this is not true. The coupling 
/33 ~ /{2miVB-L) should remain perturbative (/33 < 1). This means that for a given 

mi, there is a lower limit on vb-l- For example, if mi = 10“^ eV, and for mo ,3 = 115 
GeV, s = 0.52, it must be that vb-l > 1-8 x 10^® GeV so that < 1. It turns out that 
with this constraint, the branching ratios (BR) for the decays r —jU + 7 and /x — e + 7 
are suppressed within mSUGRA, mainly because the logarithm upon which the BR has a 
fourth power dependence, is small. We conclude that in this scenarios with 1/mi dominance 
in Eq. (16), in the case of Majorana LEV there is no signihcant LEV in lepton sector within 
mSUGRA. 

We have also examined LFV effects if the mSUGRA assumption is relaxed somewhat for 
the case of 1/mi dominance. If for example. A/ = Aq/, A^ = AqW, Ai^ = AqYi, at Mpi with 
Aq different from Aq, we found that the BR for r —/x + 7 can be in the range 10“^ — 10“® 
for relative large Aq ~ TeV. The BR for /x —> e + 7 is extremely small ( 10 “^^^) with this 

set of boundary conditions. We shall not present details of this £t, and will focus solely on 
the mSUGRA case. 

For comparison, let us also consider the case of Dirac LFV. For Mr oc 1 there is an 
unambiguous connection between the LFV effects and the neutrino mass matrix. To see 
this more explicitly, note that in this case, the seesaw formula reduces to 

-MM^ = MIMd . (19) 

Thus the product M]^Mb can be completely expressed in terms of the neutrino masses and 
mixings. Since for the GP conserving case, the flavor violating slepton mixings effects always 
involve YjY^, the LFV amplitudes can be expressed directly in terms of the neutrino mass 
matrix. Assuming again that mi m 2 m 3 , we can write down the amplitude for the 
decay r —/x + 7 to be proportional to (Am ^)23 oc while that for the decay /i —> e + 7 
is proportional to (Am^)i 2 oc csm 2 /v^- These amplitudes increase with the light neutrino 
masses, unlike the case of Majorana LFV, which scale inversely with mi. Both decays have 
branching ratios accessible in the planned experiments in the Dirac LFV alternative, as will 
be further explained in Sec. III. 


C. Enhanced LFV from Majorana Yukawa couplings 

The dominance of 1/mi in and therefore in fij is not the only possibility. We have 
found a large range of parameters where the LFV effects are significant while deviating from 
the 1/mi dominance. To see this explicitly it is more convenient to define the light neutrino 
mass matrix M.yin the gauge basis as a perturbative expansion in a small parameter e ~ 0.1 
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Ml, = mo 


( 20 ) 


( 

he™ 

de 

he™ 

1 “1“ (Z6 

1 

de 

1 

1 + be 


Here (a, b, d, e, h) are order one coefficients, ~ and the exponents (n, m) in the 

(1.1) and (1,2) entries have to be at least 1, but can be larger. This matrix provides a good 
£t to all neutrino data. If m = 1, and n = 1, 2 it will reproduce the results of the previous 
section with 1/mi dominance. But for other values of (n, m), which are also consistent 
with all neutrino oscillation data, we have found much larger LFV effects. Consider the 
case 77 , = 2, m = 4 for example. In this case the (1,1) and the (1,2) entries of Mu are not 
completely specihed by neutrino data alone. Inverting Eq. (20) we obtain for the matrix /, 

2 / (a + b)c^e^ cde^ -cde^ \ 

f = „ cde^ -d^e^ dhe^ . (21) 

d moVB-L y j 

Here we have used the parametrization Mb = Diag(ce^, e, l)mB, 3 - All parameters in / 
are hxed, except for (e, h). If (e, h) ~ 0(1), we have a situation where the off-diagonal 
entries as well as the diagonal entries of / are of the same order leading to large LFV effects. 
Contrast this with the case where m = 1 is adopted in Eq. (20), in which case we would 
have the (3,3) entry of / to be order 1, the (2,3) to be order e, the (1,3) of order e^, and the 

(1.2) entry of order e^. The choice (which is equivalent to 1/mi dominance in /) will lead 
to a suppression of the r ^ /r + 7 branching ratio by a factor 10 ^ which would put it 
in an unobservable range. The decay /i —>■ e + 7 will be suppressed even further, by a factor 

~ 10 “®, compared to the case of m = 2 . 

It is clear then that the most interesting case for lepton flavor violation is when n = 
4, m = 2 in Eq. (20). This is the choice we shall make for our numerical hts in Sec. HI. 


III. NUMERICAL RESULTS 


A. Neutrino mass fit 


We start with a basis where the charged lepton mass matrix is diagonal. For the Dirac 
neutrino mass matrix we assume that 


Mb ='jta.nldMi , ( 22 ) 

This assumption is well justihed if the model is embedded 


t/g, Tflj 


r4|. It is also possible 


where Mi = Diag(m, 

into a SUSY left-right framework as we have shown earlier papers 
that this assumption holds approximately in a wider class of models where the off-diagonal 
entries in Yi and Yu are small, although not strictly proportional. The proportionality 
constant 7 is a parameter of the model. If we assume that the Vr Dirac Yukawa coupling 
is the same as the top-quark Yukawa couplings, as happens in some versions of minimal 
S'O(IO) models, then 7=1 for large values of tan/3 ~ 60 while 7 ~ 5 — 6 for tan/3 ~ 10. 
The light neutrino Majorana mass matrix is obtained as 
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(23) 


M,.= 


7^ tan^ (3 
vb-l 


Mef-^Me . 


In our analysis, we choose values of tan /3 > 10, since smaller tan /3 values are less preferred by 
the recent Higgs mass bounds from LEP and the muon g — 2 anomaly , if interpreted 
as a SUSY effect, would prefer a moderately large tan/3. 

Following the analytic expression given in Sec. II.C, we obtain a £t for the light neutrino 
spectrum as follows. We choose tan/3 = 10 and find a good £t to the solar and atmospheric 
neutrino data with hierarchical neutrino masses if the matrix / dehned at at vb-l = 2 x 10^^ 
GeV has entries given by 


/ 


/-l.l X 10“^ 


-0.015 
V 0.29 


-0.015 0.29 \ 

0.50 -0.57 

-0.57 0.104 J 


(24) 


The resulting neutrino masses at the weak scale (using the RGE evolution 
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are: 


(mi, m 2 , ms) = (-2.7 x 10“^ 6.4 x 10■^ 8.6 x 10“^) eV . (25) 


The leptonic mixing matrix is given by: 


U 


/ 0.85 -0.52 -0.053 \ 

0.33 0.62 -0.72 

V-0.40 -0.59 -0.70 / 


(26) 


Note that the mixing angles |17e2| = 0.52 and lU/isI = 0.72 are in very good agreement 
with the solar and atmospheric oscillation data. Note also the prediction for the parameter 
jf/gsl ~ 0.05, which is in the range that can be tested in planned long baseline experiments. 

For comparison purposes we will also present LEV effects associated with Dirac Yukawa 
couplings, assuming no flavor violation in the Majorana sector. It will be desirable to have 
a neutrino mass fit for this case as well, where Mb is proportional to an identity matrix and 
all flavor mixings arise from Y^. We have found such hts for the same set of light neutrino 
spectrum as in Eq. (25) -(26). Actually, there are discrete possibilities for since the 
equations are not linear in Y^,. We have found that the predictions for LEV decays are 
essentially unchanged within these discrete possibilities. So we focus on one of them. At 
Vb-l = 9 X 10^^ GeV, we hnd that Y,y is given by 




/ 0.04 + 0.074/ 
-0.073 + 0.029/ 
V 0.025 - 0.034/ 


-0.073 + 0.029/ 
- 0.22 + 0 . 011 / 
-0.35-0.013/ 


0.025 - 0.034/ \ 
-0.35-0.013/ 
-0.24 + 0.016/ J 


(27) 


Although Y^ in Eq. (27) have complex elements, YjY^, is real, as can be easily verihed. 
Thus the light neutrino mass matrix is real (see Eq. (19)), and there is no GP violation 
in neutrino oscillations in this limit, as in the case of Majorana LEV. We shall use this £t 
for the case of Dirac LEV case in our numerical analysis. 
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mj/2[GeV] 

FIG. 1. i?(/i —> e + 7 ) vs mi /2 for three values of Aq corresponding to the case of Majorana 
LFV. The Yukawa coupling matrix / is given in Eq. (24). The horizontal solid line shows the 
current experimental limit. 

B. Predictions for lepton flavor violation 

In our calculation we make the following assumptions about the supersymmetry breaking 
sector, adopting the mSUGRA framework. We assume universal scalar masses mo, a common 
gaugino mass mi /2 and trilinear soft SUSY breaking A terms that are proportional to the 
respective Yukawa coupling matrices, with a common proportinality constant Aq. We take 
the fundamental scale at which these relations hold to be the GUT scale = 2.4 x 10^® 
GeV. We impose various experimental constraints on the SUSY spectrum. In particular, 
the Higgs boson mass limit {mh >114 GeV) along with the 6 —s + 7 constraint requires 
that mi /2 > 250 GeV. The direct limits on SUSY searches are also taken into account. 
Furthermore, radiative electroweak symmetry breaking is demanded. We choose the sign 
of /i to be positive, which is preferred by the 6 —s + 7 constraint. We do not explicitly 
impose constraints from muon g — 2, but we note that for the range of parameters we have 
chosen (positive /r and moderate to large tan/?), this anomaly can be explained via SUSY 
exchange, provided that mi /2 is limited to less than about 800 GeV (tan/? < 40). 

Additionally, we require that the lightest neutralino be the lightest supersymmetric par¬ 
ticle (LSP) in the model. Furthermore, this LSP should constitute the dark matter of the 
universe. However, for generic SUSY parameters the neutralino annihilation cross section 
is not strong enough to bring its relic abundance to the required level from dark matter 
constraints (density in dark matter should be ~30% of the critical density in order for 
the LSP to form a good dark matter candidate). This remark applies for the range of pa¬ 
rameters that we are interested in, viz., mi /2 > 300 GeV and tan/? 10. The dark matter 
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mj/2[GeV] 

FIG. 2. —> e + 7 ) vs mii 2 for two different values of tan (5. All other parameters are same 

as in Fig. 1. 


constraint can be satisfied by invoking the co-annihilation mechanism involving the LSP and 
the second lightest SUSY particle (NLSP), which is the right-handed stan in these models 
I^ . In order for this co-annihilation to be efficient, we need the lighter stan mass to be 


within 25-30 GeV of the LSP mass. We choose the nniversal scalar mass mo snch that this 
constraint on the lighter stan mass is satisfied. We allow the range 0.07 < < 0.21. 

For example, for the case of Aq = 0, mi /2 — 4.4mo will be reqnired to satisfy this constraint. 

We choose = Y^Aq, Af = YfA^, A^, = Y^Aq at the GUT scale. For the common Aq, 
we choose three different values, Aq = (0, 300, 800) GeV. We find that the branching ratio 
increases with Aq, as can be inferred from the analytic approximation, Eq. ( 8 ). 

In Fig. 1 we plot the branching ratio for /r —>■ e - 1 - 7 in the Majorana LFV case as a 
function of mi/ 2 , varying it from its lower limit of 250 GeV to 1 TeV. Here we chose Aq = 0. 
We find that the BR varies between 10“^^ and 10“®. Part of the parameter space is already 
ruled out by the non-observation of /i —^ e - 1 - 7 . This experimental lower limit on the BR is 
indicated by a thick solid line in Fig. 1. We conclude that the entire parameter space can 
be explored with the planned round of /x —> e -f 7 experiment | 2 ^ . 

In Fig. 2, we show the tan/3 dependence of the B(/x —> e -t- 7 ) for the same set of input 
parameters as in Fig. 1. Here we present results for two values of tan/3, tan/3 = 10 and 30. 
The branching ratio increases with tan/3, as anticipated from Eq. (10). As in Fig. 1, we 
demand the dark matter constraint be satisfied in Fig. 2 as well. 

In Fig. 3 we plot the branching ratio for t —>■ j for the case of Aq = 0. The solid 
line corresponds to the Majorana LFV case. For comparison, we also present in the same 
graph the r ^ /x -|- 7 BR in the case of the Dirac alternative (fit of Eq. (27)). In Fig. 4 
we plot the same quantities for a different value of Aq = 800 GeV and tan/3 = 10. We find 
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mi/2[GeV] 

FIG. 3. i?(r —> ^ + 7 ) vs 777 - 1/2 (solid line) for the same set of parameters as in Fig. 2 corre¬ 
sponding to the Majorana LFV case. The dashed line corresponds to the Dirac LFV alternative 
given in Eq. (27). 
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FIG. 4. Same as in Fig. 3, but with Aq = 800 GeV. 
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from Figs. 3 and 4 that the r —/i + 7 is typically smaller in the Majorana case compared 
to the Dirac case, although they can be comparable for some values of mi /2 where there is 
an accidental cancellation between the chargino and the neutralino diagrams in the case of 
Dirac LFV (and not for the Majorana LFV for this input choice). The BR for the case of 
Majorana LFV is < 10“® for Aq = 0, but can be as large as 2 x 10“'^ for Aq = 800 GeV. 

In Fig. 5, we show the B(/i —>6 + 7 ) as a function of mi /2 for the Majorana alternative 
(solid lines) and compare it with the Dirac case (dashed line and dot-sashed line). We have 
chosen Rq = 0 and tan /? = 10 for illustration in this case. The dashed line is drawn for the 
W shown in Eq. (27). We see that the branching ratio is large in the Dirac case. However 
if we scale down the Dirac neutrino coupling by a factor of 2.5, which does not alter the 
light neutrino spectrum (since the overall scale factor is not determined from low energy 
data), we hnd that the branching ratio (denoted by the dot-dashed line) decreases. Thus, 
depending on the overall scale factor, B{fi —> e + 7 ) can be larger in the Majorana case 
compared to the Dirac case (or vice versa). 

As Fig. 5 shows, there is some uncertainty in the BR arising from the overall scale factor. 
If we consider ratios of branching ratios, this uncertainty disappears. We dehne the ratio 


B{fi ^ e + 7 ) 
R(r ^/i + 7) 


(28) 


and study its dependence on the SUSY spectrum as well as on the Dirac/Majorana cases. 
In Figures 6 and 7, we plot r for the cases of Dirac and Majorana alternatives for Aq = 0 
(Fig. 6 ) and Aq = 800 GeV (Fig. 7) for tan/? = 10. We see from these hgures that the 
ratio depends on the supersymmetry breaking parameters Aq, mo, mi /2 etc, in a way that is 
different for these two cases. This difference shows up in almost the entire parameter space. 
We also notice that the ratio r is large for the Majorana case for large Aq. We conclude that 
with some information on the SUSY breaking parameters, measurement of r can provide 
crucial insight into the origin of the seesaw mechanism. 


IV. CONCLUSIONS 

In this paper we have studied the consequences of implementing the seesaw mechanism 
by renormalizable Yukawa couplings involving a B — L = 2 Higgs held on lepton havor 
violating processes. The Majorana masses arise as renormalizable Yukawa couplings, 
rather than bare mass terms, and therefore inhuence the renormalization group evolution of 
the soft SUSY masses of sleptons. This slepton havor violation is transmitted to the lepton 
sector and induces the rare leptonic decays p —> e + 7 and r —p + 7 . We have emphasized 
that this case is qualitatively as well as quantitatively diherent from the case where the z/jj 
Majorana masses are directly put in by hand or are induced by non-renormalizable terms. 
We have found interesting structure for the Majorana Yukawa coupling matrix / which is 
consistent with low energy neutrino oscillation data and that leads to observable lepton 
havor violations in radiative lepton decays. We have adopted a SUSY breaking scenario 
with minimal havor violation, the popular mSUGRA scheme, and have found even within 
this scheme that the radiative lepton decays are within reach of planned experiments. The 
parameter space that we have explored is consistent with acceptable abundance of neutralino 
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mj/2[GeV] 

FIG. 5. —> e + 7 ) vs mi/ 2 - The solid lines correspond to the Majorana LFV case and the 

dashed line line to the Dirac alternative. The dot-dahsed line is obtained for the Dirac LFV case, 
but with a rescaling of of Eq. (27) by a factor of 2.5. The horizontal solid line shows the current 



mj/2[GeV] 

FIG. 6 . The ratio r = e + —> /j. + -y) versus mi /2 ^0 = 0. The solid 

line correspond to the Majorana LFV case and the dashed line corresponds to the Dirac LFV 
alternative. 
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FIG. 7. Same as in Fig. 6 , but with Aq = 800 GeV. 


dark matter, as well as Higgs boson mass limit, 6 —*• s + 7 constraint, muon g — 2 constraint 
as well as the direct experimental limits on SUSY particles. 

We have compared the predictions from the Majorana Yukawa coupling related LFV with 
that related to the Dirac Yukawa coupling. The ratio r = —>• e + 7 )/i?(r —> /i+ 7 ), being 

insensitive to the overall scale of i? — L symmetry breaking, along with some information 
on the SUSY spectrum, can provide deep insight into the origin of the seesaw mechanism. 
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